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Abstract. We present a statistical method to derive the mass functions 
of open clusters using sky survey data such as the 2 Micron All Sky 
Survey (2MASS) and the Guide Star Catalogue (GSC). We have used this 
method to derive the mass functions in the stellar/substellar regime of 
three young, nearby open clusters, namely IC 348, a Orionis and Pleiades. 
The mass function in the low mass range (M< 0.50M Q ) is appreciably 
flatter than the stellar Salpeter function for all three open clusters. The 
contribution of objects below 0.5 Mq to the total mass of the cluster is 
~40% and the contribution of objects below 0.08 M to the total is ~4%. 



1. Introduction 

Recent surveys have found a significant population of low-mass stars, brown 
dwarfs and planetary-mass objects in young open clusters. Since low-mass ob- 
jects evolve little over the lifetime of the Universe, the present day mass function 
of these objects is a good representation of the Initial Mass Function (IMF). The 
mass function in this low-mass regime is however poorly known due to faintness 
of these objects and also due to uncertainty in the mass-luminosity relations. 
Low-mass objects at or below the Hydrogen Burning Mass Limit (HBML) of 
0.08 M are known to be warmer and hence more luminous when young although 
they cool rapidly and fade with age (Baraffe et al. 1998). The combination of 
youth and proximity in some open clusters make them ideal targets for searches 
of low-mass objects below the HBML particularly at infrared wavelengths. In 
the present study we have adopted a statistical approach to determine the mass 
function (dN/dM oc M — Q ) of objects in the mass range 0.5 M Q to 0.025-0.05 M 
using data of three open clusters namely IC 348, a Orionis and Pleiades. 
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2. Sample Selection and Analysis 

We have used the data from two extended sky surveys — 2MASS Second In- 
cremental Release and the latest version of GSC — on the three open clusters. 
The limiting magnitudes of 16.5, 15.5, 15 and 21 in the J, H, K and F (POSS 
II IIIa-F) passbands, respectively, enable us to probe down to 0.025 M in IC 
348 and a Orionis and O.O4M in Pleiades. Unlike most other previous studies 
which rely on confirming candidate low- mass objects through spectroscopy we 
use a statistical approach to estimate the number of low-mass objects. In this 
approach it is important to use several control fields close to each cluster to 
subtract the contribution of foreground and background objects. The nature of 
these two extended surveys permitted us to use several control fields and vari- 
able field sizes. Table 1 lists the positions and radii of the fields chosen for the 
three clusters. 



Table 1. Positions and sizes of the fields 



Fields 


RA (J2000) 
(h m s) 


Dec (J2000) 

(° ' ") 


Radius 
(arcmin) 


IC 348 


03 44 30 


+32 17 00 


20 


Control 1 


03 49 08 


+31 19 08 


20 


Control 2 


03 44 10 


+33 19 26 


20 


a Orionis 


05 38 45 


-02 36 00 


30 


Control 1 


05 58 29 


-04 29 48 


30 


Control 2 


0.5 11 00 


-00 20 00 


30 


Pleiades 


03 47 00 


+24 07 00 


90 


Control 1 


03 18 00 


+26 41 00 


90 


Control 2 


03 05 00 


+24 42 00 


90 



In the first stage of analysis we merged the 2MASS and the GSC sources by 
taking 2MASS coordinates and cross correlating them with the GSC catalogue. 
A search radius of 2"was found to be the optimum value for the cross correlation 
being small enough to reject spurious and multiple detections and large enough 
to include any minor positional uncertainties in the two catalogues. Selection 
criteria to pick out cluster members were derived by inspecting various colour 
magnitude diagrams (CMDs). The use of the F versus (F-R) CMD overcomes 
the degeneracies between the mass and distance found in the K s versus (J-K s ) 
CMD, and minimizes the overlap between reddened background stars and low- 
mass members of the cluster. This leads to a more efficient rejection criteria 
for non members which is accomplished using the empirical data from Leggett 
(1992) and the theoretical isochrones of Baraffe et al. (1998). The selection 
criteria for low-mass members, derived using the suitably distance-scaled and 
extinction corrected model isochrones of Baraffe et al. (1998), are listed in Table 
2. The sources satisfying the criteria are further statistically filtered for possible 
contaminants using the control fields. Details of the procedure are discussed in 
Tej et al. (2002). Finally, the masses of the selected candidates are determined 
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by comparing observed magnitudes with those derived from evolutionary models 
(Baraffe et al. 1998 and Chabrier et al. 2000). 



Table 2. Criteria for low-mass cluster members 



IC 348 a Orionis Pleiades 



F < [2.58(F - J) F < [2.55(F -J) F < [7.93 + 1.99(F - J) 

+ 10.14] +10.26] +056(F - J) 2 - 0.13(F - J) 2 

+0.01(F - J) 4 ] 

F — K > 4.09 F-K> 4.35 F — K > 3.42 

K > 11.20 K > 11.17 if > 11.17 

J - if > 0.90 J — if > 0.94 J - if > 0.83 

F- if > 0.20 H-K> 0.20 H-K> 0.22 



3. Results 



We first carried out a detailed study for IC 348. For this cluster, we derived 
the mass function using the solar metallicity models of Baraffe et al (1998) and 
dusty models of Chabrier et al. (2000) both of which gave similar results. The 
location of low-mass members isolated by us and the spectroscopically confirmed 
low-mass members of Luhman (1999) are in good agreement which builds the 
confidence in our selection criteria. We then used the same methodology for a 
Orionis and Pleiades. The resulting exponents of the mass functions for IC 348, 
a Orionis and Pleiades are respectively 0.7, 1.2 and 0.5 with an estimated error 
of ± 0.2. For the younger clusters IC 348 and a Orionis, there is some hint of a 
possible dip in the mass function at the position of the HBML. 



Table 3. Summary of the results 



Cluster 


Age 


Distance 


Mass Range 


a 




(Myr) 


(pc) 


(M ) 




IC 348 


5 


316 


0.5-0.035 


0.7 


a Orionis 


3 


352 


0.5-0.045 


1.2 


Pleiades 


100 


125 


0.5-0.055 


0.5 



The results of this statistical approach (Table 3) imply that though the mass 
function continues to rise above the HBML it is appreciably flatter compared 
to the Salpeter mass function. In Figure 1 we have plotted the derived mass 
function for the three clusters. Taking the canonical Salpeter exponent of 2.35 
in the mass range 1-10 M Q , the Chabrier exponent of 1.55 in the mass range 
0.5-1.0 M and the values obtained by us below 0.5 M we estimate the mass 
contribution to be about 40% for objects below 0.5 M and 4% for objects below 
the HBML of 0.08 M . Our results are consistent with previous studies (e.g. 
Bejar et al., 2001) and suggest that although low-mass objects are at least as 
numerous as their stellar counterparts their contribution to the total mass of 
the cluster is small. 
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Figure 1. The derived mass functions for the three clusters. The y-error 
bars are the yN errors involved in the counting statistics. For all three 
clusters the filled circles are points derived using the models of Baraffe et al. 
(1998) and for IC 348 the open triangles are data points derived using the 
dusty isochrones of Chabrier et al. (2000). 
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